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CAN SOLID BODY DESTRUCTION EXPLAIN ABUNDANCE DISCREPANCIES IN PLANETARY NEBULAE? 
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ABSTRACT 

In planetary nebulae, abundances of oxygen and other heavy elements derived from optical recombination 
lines are systematically higher than those derived from coUisionally excited lines. We investigate the hypothesis 
that the destruction of solid bodies may produce pockets of cool, high-metallicity gas that could explain these 
abundance discrepancies. Under the assumption of maximally efficient radiative ablation, we derive two 
fundamental constraints that the solid bodies must satisfy in order that their evaporation during the planetary 
nebula phase should generate a high enough gas phase metallicity. A local constraint implies that the bodies 
must be larger than tens of meters, while a global constraint implies that the total mass of the solid body reservoir 
must exceed a few hundredths of a solar mass. This mass greatly exceeds the mass of any population of comets 
or large debris particles expected to be found orbiting evolved low- to intermediate-mass stars. We therefore 
conclude that contemporaneous solid body destruction cannot explain the observed abundance discrepancies in 
planetary nebulae. However, similar arguments applied to the sublimation of solid bodies during the preceding 
asymptotic giant branch (AGB) phase do not lead to such a clear-cut conclusion. In this case, the required 
reservoir of volatile solids is only one ten-thousandth of a solar mass, which is comparable to the most massive 
debris disks observed around solar-type stars, implying that this mechanism may contribute to abundance 
discrepancies in at least some planetary nebulae, so long as mixing of the high metallicity gas is inefficient. 

Subject headings: atomic processes — comets: general — ISM: abundances — line: formation — planetary 
nebulae: general 



1. EVITRODUCTION 

In recent years, evidence has mounted for very large dif- 
ferences between the heavy-element abundances calculated 
in planetary nebulae (PNe) from optical recombination lines 
(ORLs) and those calculated from collisionaly excited lines 
(CELs) for C, N, O and Ne ions. Until the origin of these 
differences is completely understood, there will be a doubt 
about the accuracy of abundances derived from nebular spec- 
tra. The most comprehensive and accurate set of observational 
data has been obtained for the O^^ ion, which emits both O II 
ORLs and [O III] CELs in the visual. The ratio between the 
q2+^j_j+ abundance derived from ORLs and CEls is termed the 
abundance discrepancy factor (ADF): 



ADF = 



(0'VH+)oRL 

(02+/H+)cEL 



The presen t status of our k nowledge of ADFs is summarized 
in Robertson-Tessi & GiEnett ( 2005 ) and Liu ( 2002 2006a^. 
ADFs for the O'^"^ ion have been measured in about 90 PNe. 
Their distribution peaks at a value of 2-3, and shows a long 
tail of values exceeding 10. The ADFs tend to be larger for 
more evolved nebulae, and, in the few cases where they have 
been determined at different positions in the same nebula, they 
are larger closer to the central star. ADFs have also been found 
in HII regions, although never lar ger than about a factor of 



2 ( Garcia-Rojas & Esteban||2007 1. There is so far no fully 
convincing explanation for these ADFs, and it is not certain 
they arise from the same physical causes in the two families 
of objects. It has recently been suggested that irradiation by 
X-rays in photoionzed nebulae could produce ADFs, but the 
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available X-r ay luminosity is too small by a large factor ( Er- 
colano|2009| ). Another possibility is that the recombination 
lines arise mainly from cool, H-poor knots in the ionized gas. 
The chemical composition of those knots rules out that they 
could be made of nucleo-processes material ejected by the PN 
progenitor (Liu et al. 2000; Wesson et al. 2005). One attrac- 
tive suggestion ( Liu^2003 , 2006b,) is that they are produced 
by the evaporation o^ solid todies (asteroids, comets, planets) 
belonging to the planetary system of the PN progenitor. In 
this paper, we examine this suggestion quantitatively, by es- 
tablishing the conditions necessary for such an hypothesis to 
work. We concentrate on oxygen as the element for which 
the most reliable observational evidence exists. We first (§|2]i 
establish purely formal constraints on the size and total mass 
of the solid bodies in order that their ablation might produce 
a cool, high-metallicity gaseous phase able to reproduce the 
observed ADFs. We then (§l3]l consider in detail the physical 
processes reponsible for solidbody ablation in PNe, estimating 
the magnitude of various terms that enter into the constraints 
of §|2j In §|4]we consolidate the previous sections to derive 
quantitative lower limits on the solid body size and total reser- 
voir mass. The implications of our results are discussed in §|5] 
For convenience, a list of symbols used in the paper is given in 
the appendix. 

2. FORMAL CONSTRAINTS 

2.1. How large must the solid bodies be? 

A continuum fluid description of a collisional plasma is 
only valid for scales much larger than the mean free path for 
interactions between the particles, £. A distinct thermodynamic 
temperature can also be defined only at scales » (, since 
multiple collisions are required in order for the kinetic energies 
of the particles to reach equilibrium. Quantitative estimates of 



these length scales for nebular conditions are given in § 3.2 
Suppose that some ablation process acts on a solid body 
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of radius a so as to produce an ablative flux, Fab, of metal 
atoms, which leave the solid surface at velocity v^b- In order to 
produce a distinct high-metallicity, low-temperature phase, it 
is necessary for the ablated atoms to have a high concentration 
on scales > t. On scales smaller than {, the ablated atoms 
will stream freelyr] so their number density will be given by 
lab = F-^ToO^ Ivahr , where r is the distance from the solid body. 
If the hydrogen number density in the nebula is n, then the 
mean metallicity enhancement by number over the volume 
is given by 

zM^i"^) (1) 

\ n IP v.^\,n 

In this expression, the term F.^\,/(v.^bn) is the gas-phase metal- 
licity enhancement at the surface of the solid body, while the 
term 3a^/{^ represents the dilution over the macroscopic vol- 
ume i^, assuming ( > a. If £ < a then this term would be 
unity|2 The condition that the metal concentration must be 
sufficiently high on macroscopic scales can be expressed as 
Zabii) ^ z, where z is the "background" metallicity of the 
nebula, implying that 



a » 



> 1/2 



(2) 



This criterion establishes a minimum radius that an ablating 
solid body must satisfy in order to produce macroscopic "pock- 
ets" of metal-enriched gas of metallicity » z. In §|3]numerical 
values of this radius will be estimated for planetary nebula 
conditions. 

2.2. How much oxygen must be in the cool phase? 

In this section, we estimate what is the minimum oxygen 
mass that must be present in the cool, high metallicity gas in 
order to produce a significant ADF. This is independent of the 
origin of the high metallicity phase. 

Consider an idealized two-phase model for an ionized neb- 
ula, comprising the normal nebular gas (phase 1) plus a phase 
that is enriched in metals (phase 2). The metallicity, density, 
and temperature of each phase will be different, but we assume 
that they are in thermal pressure equilibrium. To be definite, 
we characterize the metallicity by z, the oxygen abundance by 
number relative to hydrogen. So long as z <k 1 in both phases 
and hydrogen is fully ionised, both the electron density and the 
total ion density will be roughly equal to the hydrogen density, 
n. In this approximation, thermal pressure balance implies 



niTy ^ n2T2. 



(3) 



The total oxygen mass of phase ; (/ = 1 or 2) is given by 

Mo,; = mniZiVi, (4) 

where V, is the total volume occupied by phase / and m - I6mu 
is the mass of an oxygen atom. The total luminosity (pho- 
tons s ') from phase / of a given O II optical recombination 
line is 



Li =^ ni yiZi OoRLiTi) V,-, 



(5) 



where y, is the ionization fraction of O^^ in phase / and cKorl is 
the appropriate effective recombination coefficient, which is 

^ Complications arising fromthe coupling of charged particles to the mag- 
netic field are considered in §13.21 

Equation |TJ also assumes tnat s > I, where .s is the mean separation 
between solid bodies. If this is not the case, then the right-hand side must be 
multiplied by js^ due to crowding effects. 



assumed to depend on temperature as 



-b 



(6) 



where a^. To, and - 1 are constants. 

A necessary condition for producing a significant ADF is 
that the total optical recombination line luminosity of phase 
2 should be at least equal to that of phase 1: L2 > Li. In the 
limit that most of the emission of collisionally excited metal 
lines and hydrogen recombination lines comes from phase 1 , 
this would give an ADF of at least two. Using equation Q, 
this condition can be expressed as 



ORL (T2) V2 



> 1 



or, substituting equations ([3]l, and (|6|l. 



Mo, 



> jt/min where 



f^mia — I .J, 



Numerical estimates for this minimum mass fraction, yU„ 
be provided in § [3] 



(7) 



(8) 



will 



2.3. How much mass must be in the solid body reservoir? 

The timescale for the total destruction via ablation of a solid 
body of radius a and density psoiid will be 



fab = 



n^solid 



4 3 

3 TO psolid 



4ap, 



solid 



dm^ohd/dt Tia^F-jom-^b 3Fab'«ab' 



(9) 



where m^b is the mass of each ablated atom. Thus, during the 
lifetime of the nebula, t, the total mass of ablated metals will 
be 

Mab ^ min( 1 , f /fab)M,oiid, (10) 

where Mjoiid is the initial total mass of the reservoir of solid 
bodies. 

The issue of the mixing of the ablated metals into photoion- 
ized gas is a complicated one (e.g., Stasiriska et al.,2007, ), but 
a lower limit on the required mass reservoir can be obtained 
by simply assuming that no mixing occurs, perhaps because 
the mixing timescale is greater than the nebular age. On that 
assumption, and if /o is the oxygen mass fraction of the solid 
bodies, then foM.^^, can be identified with the minimum mass 
of oxygen, Mo,2, which must be present in the cool gas phase 
in order to produce a significant ADF (see § |2.2[ ). At the 
same time, the oxygen mass, Mo,i in the background nebula is 
Mo,\ - l6z\Mi, so that the gas-phase minimum mass fraction 
condition of equation (|8]l, when combined with equation ( [T0| , 
gives the following minimum mass of the solid body reservoir: 



Msolid > 16(fab/f)Zl;"rainMi//o. 



(11) 



Note that if the solid body ablation timescale greatly exceeds 
the nebular age (fab s> f), then the mass required to be locked 
up in the solid body reservoir will be much larger th an th e 
gas-phase oxygen mass required to explain the ADF (§ 2.2 1. 



3. PHYSICAL CONSTRAINTS 

Now that we have established the formal constraints on the 
solid body size, eq. and total reservoir mass, eq. ( [TT] i, that 
must be satisfied in order for solid body destruction to explain 
the observed ADFs, we are in a position to calculate quantita- 
tive values appropriate to the physical conditions in PNe. To 
that end, we first estimate values for various parameters that 
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enter in these two equations: the minimum gas-phase oxygen 
fraction, fimm, the interaction length-scale i; and the flux and 
velocity of ablated particles, F^b and Wab- In anticipation of our 
final results, we are most interested in finding robust limits on 
these parameters such that the efi'ects of solid body destruction 
are maximised. 

In this and following sections, we frequently make estimates 
based on typical physical conditions in PNe. To show how 
the results vary with the nebular parameters, we employ di- 
mensionless scaled variables as follows: ionized gas density 
n - 10'*n4 cm"^; ionized gas temperature T - 10 T4 K; cen- 
tral star ionizing luminosity = 10^^547 s central star 
bolometric luminosity 10"* L4 Lq ; radius from central star 
R - 10'^/?i7 cm; nebular expansion velocity v - 20u20 km s"'. 

3.1. Minimum gas-phase oxygen fraction, //min 

Photoionization models indicat^ that typical values of 
T2/T1 < 0.1 are achieved if zi is enriched by a factor of 
100 to 1000 with respect to zi, which is assumed to be roughly 
at the solar value. A lower limit for {/i/{/2 can be estimated 
under the assumption that the abundance of O^ is negligible in 
both zonesFlso that the oxygen ion balance is solely between 
O^^ and O . From equation (jsjl, it can be seen that the small- 
est values of fi^in will be obtamed when there is little O"^^ in 
the cool phase 2, so that j/2 ~ 1 . To produce an appreciable 
ADF, it is necessary for the coUisional O^^ lines to come from 
phase 1, so that yi cannot be exceedingly low. We therefore 
take a conservative lower limit oi yilyi > 0.1, which, when 
combined with an estimated Ta/Ti ~ 0.1, gives a lower limit 
to the oxygen mass fraction of > 0.001. 

3.2. Particle interaction scale, { 

For the charged particles in a photoionized plasma, the inter- 
particle coupling is dominated by the accumulation of many 
small-angle deviations, mediated by the long range Coulomb 
force. These electrostatic interactions are efi'ectively shielded 
at separations larger than the Debye length, giving an effec- 
tive collisional mean free path of - 6.1 x 10^r2/«4 cm, 
where the weak dependence on physical conditions of the 
Coulomb logarithm has been suppressed. In addition, in the 
presence of magnetic fields, charged particles are strongly 
coupled to the field lines on the scale of the gyroradius: 
= 1.25 X 10^(j6Aab/n4)'^^ cm, where p is the ratio of thermal 
to magnetic pressure, which is likely to be ~ 1 but with strong 
local variations. If the magnetic fields are turbulent, then there 
may be indirect particle-particle coupling on scales ~ i]^. 

3.3. Particle ablation rate, F.^\, 

The mechanisms for the ablation of volatiles from solid 
bodies in a planetary nebula may include thermal sublimation 
(evaporation), sputtering by particle collisions, and photosput- 
tering (photodesorption) by ultraviolet photons. To be definite, 
we will consider two canonical cases for the nature of the ab- 
lating particles: (1) the case of water ice, where a relatively 
small amount of energy is required to br eak molecular bond s 
and liberate particles from the surface ( jWestley et al.|I995 I, 
and (2) the case of amorphous silicates, where much more 
energetic events are require d in order to detach particles from 
the soUd ( [Field et al.|1997 ). These will be referred to as the 
volatile case and refractory case, respectively. 



3.3.1. EUV photosputtering and the maximal ablation rate 

Photosputtering, also termed photodesorption in the case 
of volatiles, is the ejection of heavy particles (e.g., ions or 
molecules) due to the interaction of individual photons with a 
solid surface. The flux of ablated metals, Fab, will be equal to 
the incident flux of photons in the relevant spectral range multi- 
plied by 2absi'ab7 whcrc 2ab.s IS the fraction of photons that are 
absorbed by the surface and Fab is the yield. Solid bodies that 
are much larger than the photon wavelength generally have a 
low albedo, so that gabs ~ 1, but the yield is a strong function 
of the incident photon energy, being essentially zero for photon 
energies below a threshold that is of order the solid binding 
energy. The photosputtering process is similar in many ways to 
the photoelectric effect, which has been much more extensively 
studied, with the difference that it is ions instead of electrons 
that are ejected in photosputtering. This difference is crucial 
because the electrons strongly couple to the electromagnetic 
field of the photons, leading to a photoelectric yield of order 
unity for photons with energies above the electron binding 
energy of ~ 10 eV. The photosputtering yield is much more un- 
certain since no relevant laboratory or theoretical studies exist. 
It is likely to be significantly smaller than unity, at least in the 
case of refractory material, since the relatively poor coupling 
between solid ions and the photons makes it difficult for them 
to be ejected. However, a firm upper limit on the ablation rate 
can be obtained by setting Kab = L which may be a realistic 
value for volatile materials (where photochemical reactions 
can enhance the coupling) but is probably an overestimate for 
refractory solids. 

For the hot central stars of PNe, the majority of the emitted 
photons are above the hydrogen-ionizing limit of 13.6 eV, so 
that the maximal ablated flux is simply equal to the ionizing 
stellar flux: 

< = [s,l'\nR^)e-'' = 7.96 X lO" 547-^17 e"' cm"^ s"' 

(12) 

where r is the neutral hydrogen optical depth, which is small 
except for near the edge of an ionization-bounded nebula. 

3.3.2. Contribution of other ablation mechanisms 

A fundamental assumption of the above analysis is that the 
ablation occurs via a single-particle photon process and is 
limited by the stellar radiation flux impinging on the solid 
body. However, it is also necessary to consider other ablation 
processes, such as particle sputtering and thermal sublimation, 
as well as the contribution of the diffuse radiation field to the 
UV flux. 

Lyman a photosputtering — The diffuse flux of ionizing EUV 
radiation is likely to be only a smal l fraction of the direct 
stellar flux (| L6pez-Martm et al.|2001 1 and will not contribute 
significantly to the total photosputtering rate. However, non- 
ionizing far ultraviolet (FUV) radiation can also contribute 
to sputtering, especially for volatile materials such as ices, 
and this may be dominated by the diffuse component. As an 
example, we calculate the diffuse flux of the Lyman a line at 
912 A, with a photon energy of 10.2 eV. 

At typical nebular temperatures, recombination greatly ex- 
ceeds collisional excitation as a source for Lyo-, whereas, for 
Galactic metallicity, the principal sink is dust absorption rather 



' Calculations using two different photoionization codes jperland et al.| 
|1998[ |Stasinska 2005 1 give consistent results. 

reasonable assumption for the hard ionizing spectra typical of PNe. 



than escape in the line wing s (Hum mer & Kunasz|l980 Hen 
ney & Arthur|1998 Henney|2007| ). Assuming a balance be 
tween sources and sinks, the diffuse Lya flux on an opaque 
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obstacle is found to be 



Fa - nas/4o-d. 



(13) 



where is the hydrogen recombination coefficient to excited 
levels, and cr^ - 10 cm^ H"' is the dust absorption cross 
section per hydrogen nucleon in the FUV. The area that the 
solid body presents to the diffuse radiation is 4na^, as opposed 
to na^ presented to the stellar radiation, so that the ratio of the 
rate of inicidence of diffuse to stellar photons is 4F„/Ff This 
is inversely proportional to the ionization parameter, defined 
as U - Ft/cn. The value of U will vary with radius inside 
the nebula, but for an ionization-bounded, uniform density 
nebula, the mean value will be f/ = 6.2 x 10"^ {S at IR\iT4yi^ . 
Therefore, we find 



AFalF, =0.0085(r4cr_2it/)"' = 1.35/?[f ^^^'^^r^ '^V; 



-21' 

(14) 

from which it can be seen that AFa typically exceeds F», but 
by less than a factor of ten. However, in order to translate from 
Fa to an ablated particle flux F^\, it is necessary to multiply by 
the yield, which in the case of FUV photons is well-determined 
experim entally and found to be Fab - 10""' to 10"^ for water ice 
I Westley et al. 1995 ). This is much lower than the Fab - 1 that 
was assumed when calculating the maximal photosputtering 
fluxFlso that inclusion of sputtering by diffuse Lyo' would not 
significantly increase that estimate. For refractory solids, the 
FUV yields would be expected to be much lower still. 

Sputtering from pariicle impact — In general, the gas tempera- 
ture in an ionized nebula is a factor of a few lower than the 
color temperature of the radiation field. As a result, the sputter- 
ing yield from particle impacts is expected to be much lower 
than the EUV photosputtering yield. Calculations of ion sput- 
tering yields from silicates and carbon show that the yield only 
becomes appreci able for incident particle energies > 50 eV 
( [Field et al.|1997 ). In principle, this might be compensated by 
a very high impinging particle flux at the solid surface. How- 
ever, taking thermal ion speeds of ~ 10 km s ' and densities 
of ~ 10"* cm"-', it is straightforward to show that the impinging 
particle flux is generally less than the EUV photon flux of 
equation ( 12i. Therefore, particle sputtering is likely to be 



unimportant. 

Thermal sublimation — Unlike the Other processes considered 
so far, sublimation is a collective thermodynamic process that 
depends on the temperature of the solid body rather than di- 
rectly on the properties of the impinging radiation. The ra- 
diative equilibrium temperature of a solid body at a distance 
10'^/?i7 cm from a star of bolometric luminosity 10'*L4 is 



solid 



48.1 l\'*R-''Hq,UQ. 



(15) 



where Qabs and Qem are respectively the absorption and emis- 
sion efficiencies of th e solid body with respect to a black body. 
The results of § |4.I| rule out bodies smaller than meter-size, 
in which case a » A and one has Qabs ~ Gem ~ 1 • The sub- 
limation flux from the s urface of water ice is su ch a strong 
function of temperature ( [Grigorieva et al. 2007 1 that it may 
be considered a step function at Tsub - l50 K. From equa- 
tion ( |I5| l, this may be transformed into a radius in the neb- 
ula within which sublimation will dominate photosputtering: 

' It is entirely possible that the sputtering yield for ionizing EUV photons is 
also <K 1 , in which case the total ablation flux would have to be revised down- 
wards from its maximal value. However, we know of no relevant experimental 
studies that address this question. 



Rn < O.I(150/rsub)^i4^^ The bolometric luminosity of the 
central star will initially be rather high (L4 ~ 1), falling to 
L4 ~ 0.01 once envelope hydrogen burning ceases (Bloecker 
|I995| l, with the higher-luminosity phase lasting for most of 
the lifetime of the nebula if the central star mass < 0.6 Mq, 
but only a fraction of this time for higher mass central stars. 
Therefore, thermal sublimation of ices can be important lo- 
cally during the early, compact phases of PN evolution (the 
first few hundred years, assuming an expansion velocity of 
20 km s"'), but proceeds at a rate much lower than maximally 
efficient photosputtering over the vast majority of the nebular 
lifetime and volume. Refractory materials, such as silicates, 
have a sublimation temperature that is roughly 10 times higher 



( Guhathakurta & Draine|[T989| e.g.,) and so sublimation of 
these is never important beyond a few AU from the star. 

In summary, none of these ablation mechanisms is likely to 
significantly reduce the ablation lifetime from the value cal- 
culated assuming maximally efficient photosputtering (equa- 
tion ([T2|). 



3.4. Particle ablation velocity, Uab 

An upper limit to the velocity of the ablated metal ion can be 
found by assuming that its kinetic energy, E^\„ is a significant 
fraction of the incident photon energy, giving E.^\, < 30 eV. 
The velocity is then u-^b = (2-£'ab/'«ab)'^^- 

4. RESULTS 

4.1. Minimum size of ablating solid bodies 

Substituting the value for the maximal ablation rate into 
equation Q, one finds a dependence on n/f», which is pro- 
portional to the inverse of the ionizatio n pa rameter, calculated 
in § 3.3.2| Hence, using the result of § 3.4 the size constraint 
becomes 



ab 



3Fabf/ 



1/2 



= 0.0012 



z_-,HRn * 



(16) 



where Aj^ = Aab/16 and z_3 = z/10"^ (f or reference, the solar 
oxygen abundance is z_3 = 0.49 + 0.06, Asplund et al.|2009] l. 
Note that this result is rather insensitive to changes in the 
nebular parameters. 

From § 3.2 a representative value of the interaction scale is 
{ - 10^ cm, so it can be seen from equation ( 16 1 that only the 
ablation of bodies larger than ~ 20 m can produce cool, metal- 
rich pockets of ionized gas. In particular, this local constraint 
rules out dust grains, but allows larger bodies such as comets 
or planets 

4.2. Minimum mass of solid body reservoir 



Turning now to the global constraint of §§ 2.2 and 2.3 the 
nebular mass can be estimated as 



M ^ -nR'nmn ^ 0.03(5 47 7^4/? 17) '^^^( 



(17) 



In principle, there is a second regime of small sizes that can also give 
ZabiO > zo due to a very large number of evaporating bodies being present 
within the volume (see footnote Wj. However, it can be shown that this 
condition is a/{ < 3 X 10""*/S47/ljj /(psoiid^iy), where / is the ratio of the 
total mass in solid bodies to the total nebular ionized gas mass and Psolid is 
the density of each solid body (in g cm"'). It can be easily seen that this size 
regime is unphysically small for any reasonable combination of parameters. 
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so that equation ([8| for the minimum mass of ablated metals 
becomes 



> 9.0 X \0-\^^^l^^{S4lTAR\^y'^ Mg 



(18) 



Using the an ablation rate of Fab times the maximal value of 
equation ( 12 1, we find an ablation lifetime of 



fab = 1.7X 10** 



2 \ 



'fl5Psolid^l7 
.^lei'ab'S 47 ) 



years, 



(19) 



whereas the nebular age can be expressed in terms of its radius 
and a mean expansion velocity, which we normalize by a 
typical PN value of 20 km s ' : 



20m{Rnlv2ii) years. 



(20) 



Hence, the minimum mass required in the solid body reservoir 
is given by 



solid 



(21) 



5. DISCUSSION 

5.1. Ongoing destruction of solids during the planetary 
nebula phase 

The constraints derived in the previous section for destruc- 
tion of solids during the planetary nebula phase are very severe. 
Solid bodies must be meter-sized or larger, ruling out classical 
interstellar or circumstellar dust grains. Furthermore, the total 
mass of solid bodies must exceed a few hundredths of a solar 
mass for kilometer-sized bodies, scaling linearly with the solid 
body radius. These limits apply to volatile materials, such as 
water ice, where the EUV photosputtering yield, Y^b, may ap- 
proach unity. Refractory materials, such as silicates, are likely 
to have a lower yield, which would require an even higher 
reservoir mass. Any dilution of the high metallicity gas by 
mixing with the background nebula over the lifetime of the PN 
would also only serve to tighten the mass constraint. Potential 
populations of solid bodies that might be found around the 
central star of a PN include planetary systems, comets, and 
debris disks. We will now examine each of these in turn. 

Planets can be ruled out as significant sources of gas phase 
metals on several grounds. First, they are expected to be 
found rather close to the central star, < 100 AU, which is 
much smaller than typical PN sizes, which are > 10'* AU. 
Second, the most massive planets are likely to be gas giants, 
which, although they may be evaporated during the PN phase 
( [Villaver & Livio 20 07 ), are unlikely to contribute to oxygen 
enrichment. Rocky planets of size 1000-10,000 km could 
potentially provide the metals, but equation 1 21 1 with a^ - Iff* 
implies a minimum reservoir mass of ~ lODD Mq, which is 
obviously out of the question. 

Comets in the Oort cloud around our sun have semi-major 
axes th at typically lie in t he range 1,600-100,000 AU, or 0.005- 
0.5 pc ( |MorbideUi|2005[ l, which is similar to the sizes of plan- 
etary nebulae. The individual comets are mostly kilometer- 
sized, thus easily satisfying the local constraint of § 4.1 How- 
ever, the total number of comets in our solar systeni has been 
variously estimated as 10"-10'^ (Neslusan 2007), with typical 
individual masses of 4 x 10'^ g ( ^Heisler^l990) , giving a total 
mass of 10 ''-10"^ Mq for the comet reservoir. This is many 
times lower than the estimated minimum mass of ~ 0.03 Mq 



Debris disks around main sequence stars are detected via 
their FIR dust emission, but this dust is required to be con- 
stantly replenished by collisions among an unseen popula- 
tion of planetesimals (Wyatt 2008). The total mass of the 
debris disk is initially ~ 10^^-10 '* Mq, and this is likely to 
have fallen by a factor a few by the time t hat the star reaches 
the AGB and PN phases of its evolution ( Lohne et al.|2008 1. 
Taking the minimum-sized body (~ 10 m) that can satisfy 
equation ( 16i, then equation (21 1 requires a reservoir mass 



of ~ 10 Mq. This already exceeds the total planetesimal 
mass of the most massive debris disks, but the situation is 
even worse since models suggest that the debris disk mass is 
dominated by larger bodies, with the mass of rn eter-class rocks 
being as low as 10 '° Mq ( |Lohne et al.|2008| Fig. 5). Since 
the required reservoir mass scales linearly with the radius of 
the body, this would be ~ 0.1 Mq for kilometer sized bodies, 
which is again many times too large. 

5.2. Comet evaporation during prior stellar evolutionary 
phases 

The preceding sections have concentrated on the destruction 
of solid bodies during the relatively short-lived planetary neb- 
ula phase, but if a much longer span of time were available 
for the evaporation process, then a less massive solid body 
reservoir would be sufficient. The majority of the physical 



processes considered in § 3.3 are inoperative during earlier 
stages of stellar evolution, smce they depend directly or indi- 
rectly on the existence of ionizing photons. However, thermal 



sublimation (§ |3.3.2.0 i is an exception to this since it depends 
only on the bolometric luminosit y of the central star, whic h is 



generally high (lO^-lO'* Lq, e.g.," |Marigo & Girardi|2007 ]l for 
approximately 10^ years while the star ascends the thermally 
pulsing asymptotic giant branch (AGB), immediately prior to 
the formation of the planetary nebula. Indeed, observations of 
water vapor emissio n in the circ umstellar outflows of carbon- 
rich AGB stars (Melnick et al.|,200 1 ) have been interpreted 
( Ford & Neufeld|2001| i as evidence for the sublimation during 
this evolutionary stage of a population of comets lying at a 
distance from the star of ^ 100-200 AU, although current 
observations are insufficie nt to discriminate between th is and 
competing mechanisms ( Gonzalez-Alfonso et al.|2007 l. 

The high-luminosity AGB phase typically l asts 100 to 1000 
times longer than the pl anetary nebula phase ( Bloecker|1995[ 
|Marigo & Girardi|2007| l, so that the estimated minimum mass 
of the solid body reservoir (eq. ( [TT| i) is reduced accordingly, to 
~ 10"^ Mq, which is now comparable with estimates for the 
total mass of comets orbiting the sun (see above). A similar 
result was obtained from the more detailed comet sublimation 
models of Ford & Neufeld ( 2001| l. However, even if mixing 
can be avoided over such a long timespan, it is not possible 
that comets evaporated during the entire AGB phase can ex- 
plain the abundance discrepancies in planetary nebulae. This is 
because the AGB stellar wind speed (typically 5 to 25 km s"'. 



imphed by equation (21 1. 



Loup et aL]|1993 i is a significant fraction of the subsequent 
expansion speed of the planetary nebula shell (20-30 km s 
Weinberger 1989; Huggins et al. 2005), while the observed wa- 
ter vapor emission line has a similar width to CO lines from the 
same AGB star ( Gonzalez- Alfonso et al.|2007 ), implying that 
newly evaporated metals are efficiently entrained by the wind. 
Therefore, a planetary nebula of age f years will incorporate 
only those metals that were evaporated within a period of at 
most ^ 6f years before the onset of the planetary nebula phase. 
Any metals that were evaporated earlier than this, during the 
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previo us > 10^ years of AGB evolution Marigo & Girardi 
( 2007 1, will lie in the low surface brightnessl)uter halo of the 
planetary nebula ( |Balick et al.|1992t [Waver et al.|2002a|b 



and so cannot contribute to observed abundance discrepancies. 
Since the bolometric luminosity of the star towards the end 
of the AGB phase (~ 10"^ Lq) is very similar to that at the 
b eginning of the planetary nebula phase, the estimates given in 
§ [3.3.2.0 for the sublimation contribution would be increased 
by at most a factor of a few. 

To obtain a better quantitative estimate of this effect, we 
use the modeling results of For d & Neufeld (2001 ), who find 
an average water outflow rate of ^ lO^^'MsoUd yr"' during the 
late AGB phase, where Msoiy is the total ice mass in suitably 
situated bodies]^ The reciprocal of this quantity (10^ years) is 
equivalent to the ablation timescale, fa[,> of equation ( 19 1. Us- 
ing equation ( [20| for the nebular age, together with the above 
estimate that ao times longer AGB mass-loss timespan can 
contribute material to the nebula, gives an effective timespan 
for sublimation of f ^ 10"* years. From equation (18i, the 
minimum evaporated oxygen mass for significant abundance 
discrepancies is ~ 10"^ Mq, thus the required comet reservoir 
mass is f^b/f = 100 times larger, or ~ 10""^ Mq ^33 M^, for 
comets with astrocentric orbital radii in the narrow range 200- 
300 AU. Closer-in comets are completely vaporised during the 
early AGB phase, whereas farther-out comets never attain high 
enough temperatures to undergo significant sublimation. The 
estimated mass of our Solar System's entire present-day Oort 
cloud ( Neslusan 2007 ) is 10-100 times smal ler than this, while 
the mass of the "extended scattered disk" ( [Morbidelli|2005[ ) 
that occupies the required range of orb ital radii is much smaller 
still (~ 0.01 Mffi, [Gomes et al.|2008] ). However, the required 
comet reservoir mass is not rnuch larger than the 5-10 
that is requirecp^ in order for comet evaporation to explain 
the H2O 1 iQ- lpi emission line from the carbon-rich AGB star 
IRC+10216 ( [Ford & Neufeld|2001| [Gonzalez- Alfonso et al.[ 
[2007] l. 

Thus, evaporation of the same icy bodies could potentially 
explain both the presence of water vapor in a carbon-rich AGB 
stars and the abundance discrepancies that are observed in 
the subsequent planetary nebula phase, but only if the total 
volatile mass were orders of magnitude large r than any comet 
population observed in our Solar System. [Ford & Neufeld 
attempt to resolve this discrepancy by positing a much larger 
(> 10 Me) initial mass for the Kuiper belt, which was subse- 
quently ejected to larger heliocentric radii, where it remains 
unobserved. However, more recent deep surveys of Trans- 
Neptunian objects have fai led to find any evidence fo r such a 
population beyond 60 AU ( Fuentes & Holman|2008 1. On the 
other hand, it is possible that the low mass of our present-day 
Kuiper belt is due to a delayed planetesimal-driven orbital 
migration of the giant planets (Tsiganis et al. 2005), which 
destabilised the outer planetesimal disk, transferring many 
bodies into the in ner solar system, w here they collided with 
terrestrial planets ( Gomes et al.'2005), and others into the Oort 
cloud (3000-30,000 AU,[Brasser 2008). Comparison with the 
correlation bet ween age an d mid-infrared excess in extrasolar 
debris disks (B ooth et al.|2009 ) suggests that such disruptive 
scattering events are rare in other stellar systems and that many 
debris disks may retain a high mass (> 10 M^) through to the 
AGB phase. It should be noted, however, that only ^ 16% of A- 

' The peak rate during thermal pulses is ~ 10 times higher, but the integrated 
contribution of these short bursts is small. 

Assuming the same lO'' year ablation timescale as assumed above. 



K main-seq uence stars show detectable debris disks ( Trilhng 
et al.|2008] ), which renders problematic the application of this 
mechanism to all planetary nebulae. 

It is also worth pointing out that any mixing of the evapo- 
rated comet material with the AGB wind or with the photoion- 
ized nebula has the potential to decrease the efficiency of this 
mechanism. Mixing due to thermal diffusion can be calculated 
for both planetary nebula and AGB wind conditions by suppos- 
ing that a comet of size ^10 km is instantaneously evaporated 
and that the resultant high-metalicity cloud quickly reaches 
pressure and ionization balance with its surroundings. The 
cloud temperature is assumed to be 100 K, which gives a cloud 
size of ~ 10" cm in both cases since the PN thermal pressure 
{P/k ~ 10^ cnT^ K) is very similar to the ram pressure in the 
AGB wind (assuming M - 10"^ Mq yr"', V = 15 kms"', 
R = 200 A U). Using mutual diffusion coefficients calculated 
as in Oey ( 2003 ), one finds a thermal diffusion time for the 
cloud of ~ 10"' years for the ionized planetary nebula case, but 
only ~ 10^ years for the case of the neutral/molecular AGB 
wind. Diffusion is much more efficient in the AGB wind be- 
cause the collisional mean free path is much longer due to the 
lack of long-range Coulomb interactions. This timescale is of 
the same order as the sublimation timescale calculated above, 
suggesting that diffusive mixing may be important in suppress- 
ing metallicity inhomogeneities. A more detailed study of this 
point, including turbulent mixing, will be the subject of future 
work. 

6. CONCLUSIONS 

We have shown that solid body destruction during the plan- 
etary nebula phase itself cannot explain observed abundance 
discrepancies. Although sufficiently massive populations of 
solid bodies may be present around planetary nebula central 
stars, the calculated sputtering rates are too low to ablate more 
than a small fraction of this mass during the lifetime of the 
planetary nebula. On the other hand, sublimation of volatile 
bodies during the final stages of the preceding asymptotic giant 
branch phase might possibly provide enough high-metallicity 
gas-phase material to explain the abundance discrepancies, but 
only if mixing is inefficient and in systems that possess a much 
more massive population of comets than is found in our Solar 
System. The inferred mass of solid bodies is similar to that 
required in order for comet evaporation to explain observations 
of water vapor in carbon-rich AGB stars, and is at the high 
end of the observed mass distribution of debris disks around 
solar-type stars. 
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APPENDIX 
LIST OF SYMBOLS 



Q-oRL^ Effective recombination coefficient of O line [cm s~ ] 
as: Case B hydrogen recombination coefficient [cm^ s"'] 
ySab: Dimensionless velocity of ablated particle: va,/c 
fi: Plasma parameter: Pg-^s IPm-ng 

jXram- Minimum value of Mo,i/Mo 2 that gives ADF > 1 
Psohd- Bulk mass density of soUd body [g cm"^] 
CTj: Dust absorption cross-section [cm^ H"'] 
t: Optical depth to ionizing radiation 
a: Radius of solid body [cm] 

ADF: Abundance discrepancy factor: (0^7H+)orl / (0^7H+)cel 

Aab! Atomic mass of ablated particle: m^\,/mji 

fiab! Kinetic energy of ablated particle [erg] 

fo: Oxygen mass fraction of soUd body 

Fab- Ablative particle flux [atom cm"^ s"'] 

Fa- Flux of Lyman a photons [cm^^ s~'] 

F,: Ionizing photon flux [cm^^ s"'] 

i: Label for normal (i = i) and metal-enriched (i = 2) nebular gas phases 

f: Characteristic length scale for interparticle interactions [cm] 

L, : Emission line luminosity of phase i [photon s"'] 

L,: Bolometric stellar luminosity [erg s"'] 

niab: Mass of an ablated particle [g] 

'"soM- Mass of an individual solid body [g] 

mn: Mass of a hydrogen atom [g] 



m: Mass of atom or molecule [g] 
MabC Total mass of ablated metals [g] 
M: Total mass of the ionized nebula [g] 
^solid • Total mass of solid body reservoir [g] 
Moj'. Total oxygen mass of phase / [g] 
Mq: Solar mass [g] 

nab: Number density of ablated metal atoms [cm"'] 

n, : Hydrogen number density of phase i [cm~'] 

n: Hydrogen number density of nebula [cm"'] 

S , : Stellar ionizing photon luminosity [s" ' ] 

2abs '• Solid body absorption efficiency 

Qem: Solid body emission efficiency 

r: Distance of ablated particle from soUd body [cm] 

R: Distance from star, radius of nebula [cm] 

rab: Timescale for destruction of soUd body by ablation [s] 

n Age of nebula [years] 

T: Gas temperature of nebula [K] 

T, : Gas temperature of phase (■ [K] 

TsoM- Temperature of solid body [K] 

Tsub: Sublimation temperature [K] 

U: Ionization parameter 

Dab: Velocity of ablated particle [cm s~'] 

v: Expansion velocity of nebula [cm s"'] 
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y, : Volume of phase i [cm ^] 

yii Ionization fraction, 0^*/0, in phase i 

Fab ■ Ablative yield 

z: Oxygen abundance by number relative to hydrogen 



Zi'- Oxygen abundance of phase ! 

Zab: Oxygen abundance of ablated material (before mixing) 



